ON THE SINGULARITIES OF A FREE BOUNDARY 
THROUGH FOURIER EXPANSION 

-^ , JOHN ANDERSSON, HENRIK SHAHGHOLIAN, AND GEORG S. WEISS 

Abstract. In this paper we are concerned with singular points of solutions 
to the unstable free boundary problem 

^1 ^" = -X{i,>0} in ^1- 

The problem arises in applications such as solid combustion, composite mem- 
branes, climatology and fluid dynamics. 
^H ' It is known that solutions to the above problem may exhibit singularities — 

■^^ , that is points at which the second derivatives of the solution are unbounded — 

• 1 as well as degenerate points. This causes breakdown of by-now classical tcch- 

^S^ • niques. Here we introduce new ideas based on Fourier expansion of the non- 

^ \ linearity X{u>o}- 

The method turns out to have enough momentum to accomplish a complete 
description of the structure of the singular set in R"^ . 
A surprising fact in R^ is that although 

u{rx) 
>■ ' sups^ |w(rx)| 

Ch , can converge at singularities to each of the harmonic polynomials 
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it may not converge to any of the non-axially-symmetric harmonic polynomials 
lO ■ a ((1 + S)x'^ + (1 - <5)?/^ - 2^2) with S ^ 1/2. 

(_P Wc also prove the existence of stable singularities in R'^ . 

o' 



Contents 

1. Introduction 2 

2. Notation 4 

3. General Background 5 

4. Fourier Series Expansions of Global Solutions 8 

5. Growth of the Solution 13 

6. Existence of a True Three-dimensional Singularity 15 

7. Estimating u — Jl{u) — Zn(u) 16 



2000 Mathematics Subject Classification. Primary 35R35, Secondary 35B40, 35J60. 

Key words and phrases. Free boundary, regularity of the singular set, unique tangent cones, 
partial regularity. 

H. Shahgholian has been supported in part by the Swedish Research Council. G.S. Weiss has 
been partially supported by the Grant-in-Aid 21540211 of the Japanese Ministry of Education, 
Culture, Sports, Science and Technology. He also thanks the Knut och Alice Wallenberg foun- 
dation for a visiting appointment to KTH. Both J. Andersson and G.S. Weiss thank the Goran 
Gustafsson Foundation for visiting appointments to KTH. 

1 



J. ANDERSSON, H. SHAHGHOLIAN, AND G.S. WEISS 

8. Classification of Blow-up Limits in M.^ 

— An Unexpected Symmetrization Effect 18 

9. Unique Tangent Cones at True Three Dimensional Singularities 19 

10. Stable Cones 24 

11. Unique Tangent Cones at Unstable Codimension 2 Singularities 25 

12. Structure of the Singular Set in R^ 29 

13. Appendix 30 
References 38 



1. Introduction 

We investigate the singular points of solutions of the unstable free boundary 
problem 

(1.1) Am = -X{«>o} in Bi, 

arising in solid combustion (see the references in [15]), the composite membrane 
problem ([9], [8], [4], [19], [10], [11]), climatology ([12]) and fluid dynamics ([1]). 

The minus sign on the right-hand side drastically changes the problem from 
the well-known obstacle problem (see for example [5] , [7] and [6] ) into an unstable 
problem exhibiting non-uniqueness, bifurcations, unbounded second derivatives and 
more. Let us describe some of the known results. 

From standard elliptic regularity theory it follows that if u is a solution to (1.1) 
then u G C^'" for all a < 1. However, in contrast to the well-known obstacle 
problem Au = X{u>o}: the solutions to (1.1) are not C^'^ in general. The existence 
of non- regular solutions was first shown in [3] . 

For convenience let us denote the set of singular points by 

5*" = {x : w ^ C^-\Br{:>i.)) for any r > 0}. 

As observed in [15], {u = 0} is analytic and u £ C^'^ in a neighborhood of each 
X G {m = 0} n {Vu ^ 0}. Thus the set of singular points is contained in the set 
where both u and jViij vanish. 

We may expect that for x" € S*" the blow-up 

(1.2) hm "(^" + "")^, 

r->0 SUPs^(xO) juj 

should give us some information about the singular set. It was shown in [15] (see 
also Proposition 3.2 below) that at a singular point x° 

u(rjX + x°) 

lim —^ pr =P x), 

J^ooSUPs (^0) \u\ 

where p is a second order homogeneous harmonic polynomial. This raises several 
questions. 

(i) Does p depend on the choice of the sequence rj — >• 0? 
(ii) Does every second order homogeneous harmonic polynomial p appear as 

limit? 
(iii) Is there any partial regularity of the singular set? 
(iv) Do energy minimising singularities exist? 



Concerning uniqueness of blow-up limits it has been shown in [15] that in two 
dimensions the free boundary of the minimal solution close to points where the 
second derivative is unbounded, consists of four Lipschitz graphs meeting at right 
angles. In [2] this fact has been extended to any solution in two dimensions, proving 
also uniformity and quantitative estimates by methods closely related to those in 
the present paper. Concerning question (iv) it has been proved in [15] that the 
singularity in two dimensions is unstable in the sense that the second variation of 
the energy is negative. As to stability of higher dimensional singularities there 
is a gap in the proof of [15] which has been pointed out by Carlos Kenig-Sagun 
Chanillo-Tung To ([11]). In the present paper we will prove the following main 
results which among other things close the gap in [15] by showing that the cut-off 
dimension concerning this problem is 3, that is, there exist stable singularities in 
three dimensions. 
Main results: 

(i) Existence of a true three-dimensional singularity (Corollary 6.1). 
(ii) Axial symmetry of blow-up limits in three dimensions (Theorem 8.1). 
(iii) Unique tangent cones at true three-dimensional singularities in R^ (Theo- 
rem 9.1). 
(iv) Unique tangent cones at unstable codimension two singularities in M? (The- 
orem 11.1). 
(v) Stability of true three-dimensional singularities in M? (Theorem 10.1). 
(vi) Regularity of the singular set in three dimensions (Section 12). 

Discussion. In contrast to the analysis of singularities for minimisers or stable 
solutions, where there are many methods available, there are few results on unique 
tangent cones at unstable singularities. Even the Lojasiewicz inequality approach 
(see for example [20]J would be hard to realize in our problem due to the lack of a 
suitable local Lyapunov functional; we do have a monotonicity formula playing the 
role of a local Lyapunov functional, but as it turns out it has the wrong scaling to 
be used at the unstable singularities of "supercharacteristic growth". 

The natural approach would be to study blow-up limits in order to analyze the 
singularities. Unfortunately the blow-up sequence in (1.2) does not provide enough 
information of the solution as the nonlinearity of equation (1.1) vanishes in the 
limit. To preserve some information of the nonlinearity we will instead, in Section 
3, consider 

(1.3) !f(!i2i±2E!)_n(.,r„x«), 

where Ii(u,rj,'xP) is the projection of u{rjX. -\- :x.^) /r'^ in Bi onto the homogeneous 
harmonic second order polynomials (see Definition 3.5). 
It can be shown that if 

u{rjy: -\- yp) 



lim — =_p(x), 

j-^ooSUPb^^.(^O) |w| 



then 



, w(rjX + x") \ 

:^™ I ^3 n(M, r.j,x.') \ = Zp, 



j^oc \ r^ 
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where Zp is a solution of 

AZp = -X{p>o}- 
Next we notice that, at each singular point x*^, 

n(u,r,-,x°) ,. u(r,x + x'') 

lim = lim . 

j^oo sups^ \n{u,rj,x°)\ j^oo supb_,^(xO) |u| 

So in order to prove uniqueness ofp it is sufficient to control how n(u, r, x*^) changes 
when r varies. More precisely we would want to estimate 

n(u,r,x°) n(u,r/2,x°) 



^^■^^ supsjn(M,r,x")| supsjn(M,r/2,xO)| 

Our method of proof is based on the observation that M(rx + x*^) w TrPr + Zp^ 
in Br, where Pr is a second order harmonic polynomial of norm 1. It follows that 

n(w,r/2,xO) « Il{TrPr + Zp,^, 1/2,0) = Il{TrPr, 1/2, 0) +Il{Zp,^, 1/2,0) = TrPr + 

Il{Zp^, 1/2, 0) (cf. Section 7). Therefore it is essential to control Il{Zp^, ■) in order 
to estimate (1-4)- This control will be achieved by means of an explicit calculation 
of the Fourier coefficients of Zp^ . 

Plan of the paper. 

In Section 3 we will remind ourselves of results and definitions of [22], [15], [3] and 
[2] that are relevant to the present paper. 

In Section 4 we use techniques developed in [14] based on Fourier coefficients to 
analyze Zp. We also explicitly calculate Zp whenp = 2xz ot p = ±((a;^+y^)/2— z^). 
Using these calculations we are able to show in Section 5 that in three dimensions 
and for small r 

sup|u(x + x°)| > cr^\ log r\ 

Br 

(Corohary 5.3). 

Based on this estimate on the growth of u we prove in Section 6 existence of a 
true three-dimensional singularity. 

Section 7 provides estimates on u — n(w) — ^ii(m) which we combine in Section 
8 carefully with the above analysis of Zp to show axial symmetry of blow-up limits 
in three dimensions. This is a remarkable symmetrization effect in view of the fact 
that there are of course second order homogeneous harmonic polynomials that are 
not axially symmetric. 

In Sections 9 and 1 1 we prove — once more carefully using the information gained 
on the Fourier coefficients — uniqueness of the blow-up limits at singular points in 
R"^. Based on the asymptotics in Section 9 we are able to show in Section 10 
that true three-dimensional singularities are stable. In Section 12 we use standard 
techniques to show that in the three-dimensional case the singular set may be 
decomposed into a countable set of isolated points and a component that is locally 
contained in a C^-curve. In the Appendix we have gathered technical calculations 
which may well be considered to be the core of our paper. 

2. Notation 

Throughout this article M" will be equipped with the Euclidean inner product 
X • y and the induced norm |x| . Moreover A : B = J27i=i O'ij^ij shall denote 
the inner product of two {n, n) matrices. We will use the set Q of all orthogonal 
matrices in M". We define Brix.^) as the open n-dimensional ball of center x° , 



radius r and volume r" a;„, and i?+(x") := {x e -Br(x") : x„ > 0},i3,r(x") := 
{x e i?r(x°) : Xn < 0}. When not specified, x° is assumed to be 0. We shall often 
use abbreviations for inverse images like {w > 0} := {x € 51 : u(x) > 0} , {x„ > 
0} := {x G M" : a;„ > 0} etc. and occasionally we shall employ the decomposition 
X = {xi, . . . ,Xn) of a vector x e M" . Wc will use the n-dimcnsional Lebesgue 
measure £" and the fc-dimensional HausdorfF measure 'H*''. When considering a set 
A , XA shall stand for the characteristic function of A , while u shall typically denote 
the outward normal to a given boundary. We will use Landau's symbols as signed 
variables. For example — o(l) will mean a negative quantity that turns to zero. By 
¥2 we will denote the space of second order homogeneous harmonic polynomials 
in M". We shall also use the projection 11 onto P2 as well as the norm t of n(w), 
both defined in Definition 3.5, as well as the parametrization parameters S {v) 
and S^{v) defined in Definition 4.2. Last, we shall use for p £ P2 the Newtonian 
potential Zp, i.e. the unique solution of 

AZp = -X{p>o} ill K", 
Zp(0) = |VZp(0)|=0, 

lim , ,„ = and 

IxKoo |x|3 

U{Zp,l)=0 (cf. [2] Section 4). 

3. General Background 

In this section we will gather some results from [15] and [3], and describe some 
compactness properties of blow-ups of solutions. First we will remind ourselves of 
the monotonicity formula proved in [22] . The roots of those monotonicity formulas 
are harmonic mappings ([18], [17]) and blow-up ([16]). 

Theorem 3.1. Suppose that u is a solution of (1.1) in J7 and that Bs{x.'^) C fl . 
Then for all < p < a < S the function 

*x«W ■=r-"-^ [ (\\7uf - 2max(M,0) 

-2r-"^3 f u^ dW-^ , 

JdBr(x°) 
defined in (0, 5) , satisfies the monotonicity formula 

$;^o(a) - <^^oip) = / r-"-2 f 2(Vu-zy-2-)' dW^ dr >0 . 

Jp JaB^(xO) ^ f '' 

This energy monotonicity is important since it helps us to distinguish different 
points of the set {u{0) = |Vu(0)| = 0}. In particular we may according to the 
following Proposition define the singular set S'" as 

5" = {x e Bi : w(x) = |Vw(x)| = and lim '^^{r) = -00}. 

Proposition 3.2 (Proposition 5.1 in [15]). Let u be a solution of (1.1) in D, and 

let us consider a point x*^ e 17 n {w = 0} n {Vw = 0}. 

(i) In the case $J^o(0+) = —00, lim^-i-o '' '^ " /« n / o-i ^^ d'H'^^^ = +cx), and for 
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T(x«,r):= ri-"/^„.„.u2dH"-^ 



each limit of 
u{x'^ + rx) 



as r — > belongs to ¥2 ■ 

(a) In the case $^o(0+) G (-00, 0), 



u(x° + rx) 



is bounded in W^'^{Bi{0)), and each limit as r ^ is a homogeneous solution of 
degree 2. 



^;//$"o(o+) 



0, then 
u(x°-f- 



rx) 



m W^^^{Bi{0)) asr^Q 



Remark 3.3. In [15, Lemma 5.2] it says that case (ii) in Proposition 3.2 does 
not occur in R^. Unfortunately the authors omitted the following homogeneous 
solution of second order in M^ : 

Using polar coordinates (r, 0), let 

-^cos2(0) 



u(r, </)) = 



373 -^(2<^) 



27r/3) 
' 27r/3) 

47r/3) 
■ 47r/3) 



when (j) G (7r/3, 7r/2), 
when (j) e (— 7r/6,7r/3), 
when e (— 7r/3, — 7r/6), 
when (f) G (— 57r/6, — 7r/3), 
when (f) G (— tt, — 57r/6), 
when (j) £ (— TT, — 37r/2). 



_^cos(20)-^sin(20) 
-^ cos2(0 + 27r/3) + ^ sin(2(/) - 
-2^ cos(20 + 2^/3) - ^ sin(2(^ ■ 
-^ cos2(</. + 4^/3) + ^ sin(2(/) - 
-^cos(20 + 4V3)-^sin(2(/)- 

Then m is a second order homogeneous solution to equation (1.1). 

Let us show that up to rotations, u is the unique non-trivial second order ho- 
mogeneous solution to equation (1.1) in R^. Each cone in which u is negative has 
to have an opening of exactly 7r/2. Thus u can be negative in at most 4 different 
connected components. However if u is negative in four components then u < 
and thus Au = 0. Since m(0) = 0, w = by the strong maximum principle. If 
u is negative in only one component then Am = — 1 in a cone with opening 37r/2 
with zero boundary values on that cone, and it is easy to see that such a u is not 
homogeneous of second order. If it < in two components then, after a rotation, 
Au = —X{xy>o}- In "2-dimensional Solutions" (vi) in Section 4, we will see that 
such a solution is not homogeneous either. The only remaining possibility is that 
M < in three components, with three components of u > in between. As the 
gradient of u is continuous across the zero level set and u is symmetric in each cone 
where u has a sign it follows that the opening of the cones where u > must equal 
each other. It follows that u > in three cones of opening tt/6 where Au = — 1. 
Thus u is unique up to a rotation. 

In [3] the authors have obtained existence of solutions in two dimensions ex- 
hibiting cross-like singularities at which the second derivatives of the solution are 
unbounded (case (i) of Proposition 3.2), as well as degenerate singularities at which 
the solution decays to zero faster than any quadratic polynomial (case (iii) of Propo- 
sition 3.2): 



Theorem 3.4 (Cross-shaped singularity, Corollary 4.2 in [3]). There exists a so- 
lution u of 

Au = -X{u>o} i'n Bi C K^ 
that is not of class C^'^ . Each limit of 

u(rx) 

r(M 

as r — > coincides after rotation with the function (x\ — X2)/\\x\ — 2^211^^(031 (o))- 

Remark on the Proof: In [3] the authors show that one can construct a solution 
u to (1.1) with linir^o *o ('') ^ ~^ ^r any M > 0. Then they use [15, Lemma 5.2] 
that states that if limr^.o 'I'o (^) < in R^ then we are in case (z) of Proposition 3.2. 
As we pointed out in Remark 3.3, Lemma 5.2 in [15] is not true. The proof in [3] 
however is easily fixed: we only have to notice that all second order homogeneous 
solutions of (1.1) have fixed energy $q = toq which follows from the uniqueness 
in Remark 3.3. Thus if we choose the constant M large enough we can exclude 
the possibility that we are in case (ii) or (iii) of Proposition 3.2 and the theorem 
follows. D 

The proof of the previous theorem can be adjusted to construct other kinds of 
singular points (see Corollary 6.1). 

Definition 3.5. By n(u,r, x") we will denote the projection operator onto P2 
defined as follows: n(u, r, x°) = rp, where r e R+ and p e P2 satisfies sup^^ \p\ = 1 
as well as 



inf 



D2M(rx + x°) „2, ^_ f D'^u{rx 



-D^h 



-tD^P 



2 



'Bi(O) 

Wc will often write n(u, r) when x*^ is cither the origin or given by the context. 
At times wc will also denote t^ ~ sup^^ |n(u,r)| and pr = n(u,r)/r,-. 

The following Lemma justifies the previous Definition. 

Lemma 3.6. The following four statements hold. 

(i) For each v G W^'^{Bi) the minimizer of Definition 3.5 exists and is unique. 

Thus n : W'^^'^{Bi) x (0,s) x Bi_,, -^ P is well-defined. 

(ii) n is a linear operator. 

(iii) If h £ W^'^{Bi) is harmonic in Bi then n(/i(x), s) = n(/i(x), r) for all s, r G 

(0,1). 

(iv) For every v,w G W^'^{Bi), 

sup |n(u + w,r)\ < sup |n(w, r)\ + sup |n(i(;, r)\. 

B\ B\ B\ 

Proof. The first and second statement follow from the projection theorem with 
respect to the L^(i3i; M" )-inner product and the linear subspacc {/ £ L-^(i3i; E" ) : 
/(x) is symmetric, constant and tracc(/(x)) — 0}. 

Writing h as the sum of homogeneous harmonic polynomials hj that are orthogonal 
to each other with respect to 



{v,w) := V dijvdijW, 

we see that n(/ij) = for all j such that the degree of hj is different from 2, 
implying the third statement. 
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The last statement follows from the linearity of 11 and the triangle inequality in 
L2(Bi;R"'). D 

Next we mention that solutions to (1.1) have second derivatives in BMO. This 
has been proved in [2, Lemma 5.1] using standard facts of harmonic analysis. 

Proposition 3.7 (cf. [2, Lemma 5.1]). Let u he a solution to (1-1) in Bi such that 
supsi \u\ < M and m(0) = |Vm(0)| = 0. Then 



sup 

Bi 



— 5 n(u,r) 



< C'o for every r < — 



where the constant Co depends only on M and n. 
Furthermore, for each a < 1 

||^^-n(u,r)||^,,„(^^ <a(M,n), 

and for each p < oo 

||A_i_nKr)||^,.,(^^)<Cp(Af,n). 

4. Fourier Series Expansions of Global Solutions 

In this section we will remind ourselves of the work of L. Karp and A.S. Margulis 
[14]. In particular, Theorem 3.1 and Proposition 3.2 in [14], summarized in the next 
theorem, will be of importance to us. 

Theorem 4.1 ([14]). Let a £ L°°(R") be homogeneous of zeroth order, that is 
CF^x.) — CT(rx) for all r > 0. Assume that a has the Fourier series expansion 



(j{x.) = y^ajCTi, 



on the unit sphere, where at is a homogeneous harmonic polynomial of order i. 

Moreover assume that Am = a and that u(0) = | Vm(0)| = linix-i-oo ■'^t(x)/|x|^ = 0. 
Then 

M(x)=(?(x)log|x| + |x|20(x), 
where 

0.2 

1= — —^^2 



and 

|2 



X 



-- '^l"g(n + ^)(^-2)"^^|x|^ 

Let us explain how we are going to use Theorem 4.1 in the present paper: If u 
is a solution to equation (1.1) such that u(0) — |Vu(0)| = 0, if 

u(rjx) 

lim , = p 

i-i-oo sup^^ \u\ 



for some p £ P2 and some sequence r^ — ?► 0, and if 




then by C^'"-convergence we will have Zp{Q) — |VZp(0)| = as well as AZp = 
— X{p>o}- Also, by weak VF^'^-convergence we will have 

= / D^Zp : D^h 

Jbi 

for all /i e P2. The latter is equivalent to n(Zp, 1) =0. 
By Theorem 4.1 with a ~ ~X{p>o} "^c can write 

Zp(x) = g(x)log|x| + |x|2</,(x), 

where q = i^r^i'^i- From here on we will assume that ri = 3 and x = (x, y, z). It 
will also be convenient to parametrize the second order harmonic polynomials. We 
will assume that p = p5 ;= (l/2 + (5)a;^ + (l/2 — (5)y^ — 2;^. This can be done without 
loss of generality since there is always a rotation of the coordinate system such 
that D^y is a diagonal matrix. Rotating the coordinate system in that way, and if 
necessary renaming x, y and z wc can always make sure that p is up to a scaling 
factor of the form above or that — p = (l/2 + (5)a;^ + (l/2 — (5)y^ — z^. The latter case 
can be handled similarly. We would want to calculate 02- To that end we choose 
the polynomials "ix^ — |xp, 3y^ — |xp and 3z^ — |xp spanning the axisymmetric 
second order harmonic polynomials in M^. That choice is somewhat arbitrary, but 
we contend that choosing different polynomials would not facilitate substantially 
anything that follows. It follows that 

02 = C((3A,((5) - A{5))x^ + {iAy{b) - A{b))v' + (3A,(5) - A{8))z'), 

where C has been chosen such that ||cr||L2(asi) = 1- 

Using spherical coordinates a; = r sin(^) cos(0), y = r sin(^) sin((/)) and z = 
rcos(0), the characteristic function X{p,>o} = X{e>arccot(V i/2+^cos(20) )}- ^^e coef- 
ficients satisfy 

JdBi Jo Jarccot(-y/l/2+i5 cos(20)) 

M^)^ f -X{ps>o}y' ^ -8 f f , Sin3(0)sin2(</>)d0#, 

JdBi Jo Jarccot(-^l/2+i5 cos(20)) 

A,{6)^ -X{p,>o}Z^ ^ -8 / sin(0)cos2(0)d0d0 



dSi Jo Jarccot(-y/l72+5cos(Z 



and 



^W = / -X{p,>o} = -8 / / , siu' 

JdBi Jo JaIccot(^Jl/2+&cos(2<t>)) 

Next we notice that with 

P,^_ 2 log 2 

5||3a;2 - l||^2(g5 x 
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(4.1) 

n(Zp,, 1/2) = -Ko((3A,((5) - A{5))x^ + {iAy{5) - A{5))y^ + {iA,{5) - A{5))z^), 

since Ii{(Ji) = for i 7^ 2. Calculating A^., Ay, A^ and A we may estimate the 
rotation of n(u,r) as follows: If u(rx)/r'^ — Il{u,r) w Zp^ then 

n(w,r)-n(u,r/2)«-n(Zp,,l/2). 

Later on this will be our main tool to analyze singular points. 

For convenience we will later also use the alternative representation ps = (1 — 
(5)x^ + 6y^ — z^ leading to the coefRcients 



BAS) = 


^7r/2 i-n/2 

-8 / / sin^(6')cos((/))^d6ld(/>, 

Jo Jarccot(^(l-5)cos2(<^)+5sin2(,/,))) 


ByiS) = 


^7r/2 ^7r/2 

-8 / / sm^(e)sm((l)?ddd(t> 



and 

I-7V/2 .7r/2 

B(d^) = -8 / / sin{e)ded(j). 

Jo Jarccot(-y/(l-i5)cos2(0)+5sin2(0))) 

That is, i?a;(l/2 — (5) = Ax{S) etc. It will be convenient to define the parameter 5 
for polynomials and solutions: 

Definition 4.2. For each v g VF^'^(i3i), let, if necessary after rotation, 

U{v, 1)/ sup |n(w, 1)1 =(1/2 + (5).t2 + (1/2 - S)y^ - z^ 

or - (1/2 + 6)x^ - (1/2 - 6)y^ + z^. 

We note that S is unique and define S^{v) := (5. 
Moreover, let 

n(w, 1)/| supn(u, 1)1 =(1 - S)x^ + dy^ - z^ 

Si 

or - (1 - d)x'^ - Sy"^ + z^. 

We note that S is unique and define S^{v) := S. It is important to note that 
supsi |n(w,l)|/supB^ |n(w,l)| = 1 implies S^{v) = 5>0. 
We will also use 6:^ := 6^{u{r-}) and (5f := 6^{u{r-)). 

In general we cannot calculate the integrals Ax,Ay,... explicitly. In some special 
cases however, when we have sufficient symmetry, we may even write down explicit 
solutions to the equation Au = ~X{p>o}- Luckily and surprisingly, as seen in 
Section 8, these special solutions are the only solutions appearing as limits of 

■u(rx + x°) 
-^ n(u,r,x ). 

1. 2-diniensional Solutions (cf. [2, Lemma 4.4]): 
Define v : (0, +00) x [0, +00) ^ R by 

v{x,z) := -4a;zlog(x^ +z^) + 2(.t^ - z^) (^ - 2 arctan ('-')') - n{x'^ + z^). 
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Moreover, let 



and define 



v(x, z), xz > 0,x ^ 0, 

wix, z) := ^ —v{—x, z), X < 0, z > 0, 
-u(x, z), a; > 0, z < 0, 



Z{x, z) 



w{x, z) — 7r(a;^ + z^) + 8xz 



8tt 
It has been shown in [2, Lemma 4.4] that 
(i) AZ^-X{-,z>o} inM3. 
(ii)Z(0) = |VJ(0)| = 0. 
(in) lim|x|^oo ^^ == 0. 
(iv) n(Z, 1) = 0. 

(v) UiZ, 1/2) = log(2)xz/^, t(Z, 1/2) = log(2)/(27r). 
(vi) Z is the unique function satisfying (i)-(iv). 

2. True 3D Solutions: 

Next we arc going to calculate Zp for p(x) — {x'^ + j/^)/2 — z^. 
Let us denote 

vi = 2p(x) log (x^ + y^) - 4z^ 

and 

3z|x| . 1 _, /\^ 

.|x| +z. 

Then Ai;i = and Awa = in ]R^\{x2 + ?;2 ^ q}. Also notice that 52Wi(x,2/, 0) = 0. 
Let 

w(x) = — (4w2(x)-wi(x)j 

in {p(x) < 0} n {z > 0}, where the coefhcients for vi and V2 are chosen such that 
the singularities cancel at a; = ^y = 0. Moreover, let 

^'W = - -^M^) 

3 + V31og(2-V3) , ^ |x| 
18 ^^ ^ 6 

in {p(x) > 0} n {z > 0}. 

Next we reflect t; at {z = 0} according to 

v{x, y, z) if z > 0, 



3z|x| 1 , > , /Ixl - z\ 



ZlM = ^ / \ -r n 

' v(x, y, —z) it Z < U. 

Last, we define Z\^ Z\ — n(Zi, 1) and Zi = —Z\. 
We have thus established the following lemma: 

Lemma 4.3. Let Z , Z\ and Z-i be as above. Then, with p(x) = (x^ + 2/^)/2 — z^, 
(i) AZ = -X{,z>o}, Z{0) = |VZ(0)| ^ 0, limx^^|Z(x)|/|x|3 = and 

n(Z, 1) = 0, n(Z, 1/2) = (log(2)/7r)a;z. 
(ii) AZi = -X{p(x)>o}, ^i(O) = |VZi(0)| = 0, 
limx^oo|^i(x)|/|x|3 = 
and U{Zi, 1) = 0,n(Zi, 1/2) = log(2)(V3/9)p(x). 
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(iii) AZ2 = -X{p(x)<o}, ^2(0) - |VZ2(0)| = 0, 
limx^oo |Z2(x)|/|x|3 = 
and n(Z2, 1) = 0,n(Z2, 1/2) = - log(2)(y3/9)p(x). 

Proof. The proof follows from simple calculation. 



D 



Remark 4.4. The fact that Il{Zp, 1/2) is a multiple of the polynomial p in the 
above three cases, natural though it may be, will be of paramount importance in 
later chapters when it comes to the question of unique tangent cones. 

The following two collections of properties of the A's and B's visualized in Figure 
1-5 arc of central importance in our paper and will be proved in the Appendix 
together with Lemma 4.6 below. 




0.0 0. 



Figure 1. 3Ay{S)- 
A{S) (upper graph) 

versus 3A^, ((5) -yl((5) 



Figure 

„ 3Ay{S)-Ai5} 
^- 3A,(5)-.1(<5) 

(upper graph) versus 

l-2i5 
1+2(5 



Theorem 4.5. For b e (0, 1/2); 

(z) i3Ay{S) - AiS))' > 0. 

(li) {2,A^{5) - A{6))' <0. 

(iii) 3A^(0) - A(0) = ■iAy{Q) - A(0) < 0. 

(iv) 3A^((5)- A((5) <0. 

(v) 3^,(1/2) - A{l/2) = 0. 

(vi) 3(A;' ~ A^') + 2,5(3A;' + 3A^' - 2A") + 2(3A; + iA'^ - 2A') > 0. 

(vii) 

26 



SAy{d) - A{6) 1 
3A46) - A{6) ^ 1 + 2S' 



(via) 3A^(0)--^'(0) 
M3A'y{0)-A'iO)^ 



3^3' 



47r 
3V3' 



Lemma 4.6. For the positive universal constant rjo defined in (13.6) and every 
5 €(0,1/2), 

sup |Cp5 + n(Zp,, 1/2)1 >C + 77o, 

Bi 

for every sufficiently large constant C > 0. 
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0.25 0.00 0.05 



0.10 0.15 0.20 0.25 



Figure 3. 1/Co vs. l/q 




Figure 4. Cy/Co 



0.10 0.15 0.20 0.25 



Figure 5. C^{5)/5 



Theorem 4.7. For small 5 > Q: 

(i) By {5) = By{0) - Cy{S)S, where < h < Cy{S)/\\ogS\ < K2 < +c 

(it) B{S) = B{0) - CoiS)S, where lim^^o Co{S) = +00. 

(m) bJ{5) = B^iQ) + o{5) as5^Q. 



Cy(S) 



(iv) lima^o c^ 



1. 



5. Growth of the Solution 

Since Z, Zi, Zi have growth |xp| log(|x|)| away from the origm and we expect u 
(up to some harmonic part) to be close to Zp, u should share the same growth. We 
will prove this in the next lemma. 

Before we state the lemma let us point out a simple fact that will be used 
frequently in what follows. By Proposition 3.7 wc know that if the origin is a 
singular point then u{r-K)/r'^ — n(u, r) is uniformly bounded by a constant Cq 
depending on n and ||ii||ioo(5^). This implies that when u(rx.)/r'^ is large, say 
max(sup^ |n(u, r)|, sup^ |u(rx)/r^|) > 2Co, then 



(5.1) 



- sup |n(u,r)| < sup |u(7'x)/r^| < 2 sup |n(u,r) 

2 Si Bi Bi 



So controlling the size of n(u, r) is equivalent to controlling the size of u(rx)/r^ at 
singular points. 

In the two-dimensional case the following lemma has been proved in [2, Lemma 

5.5]. 
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Lemma 5.1. Let n = 3 and let u be a solution to (1.1) in Bi such that sup^^ l"*^! ^ 
Af and u{0) = |Vw(0)| = 0. Then there exist po > and tq > such that if 

(5.2) sup|n(M,r)| > — 

-Bi Po 



for an r < j'q then 



sup \U{u, r/2)\ > sup |n(u, r)\ + ?7o/2, 



where rjo is the positive constant in Lemma 4-6 (see (13.6) in the Appendix). 

Proof. If the Lemma is not true, then there exists a sequence u^ of solutions to 
(1.1) and rj — > such that 

sup |n(w', rj)| > j and sup |n(u-', rj/2)| < sup |n(M-', rj)| + ?7o/2. 

B\ B\ B\ 

Using Proposition 3.7, and passing if necessary to a subsequence, 

^.(x) := !f!^ - ^(^.^r,) ^ v in C,t':(M^) H <f (M^). 

3 

We also have IV{u^,rj/2) = IV(u^ ,rj) + n(w^l/2). The limit v satisfies Au = 
~X{p>o}! where — passing if necessary to another subsequence — 

p^ lim p.j forpj :=n(u^rj)/||n(u^rj)||ioo(B^). 

It follows that n(ti, 1/2) = n(Zp, 1/2) where Zp is the unique solution to 

A^p = -X{p>o} inK^, 

Zp(0) = |VZp(0)| = n(Zp, 1) = and 
lim, , I^^WI - 

Consequently limj_j.oo(n(u-', rj/2) — \i{u^ ,rj)) ~ Il{Zp, 1/2), and 

(5.3) sup|n(u-'>j/2)| = sup |n(w^rj)+n(Zp^, 1/2)1 +o(l) as j ^ cx3. 

Finally wc apply Lemma 4.6 and obtain the statement of the lemma. D 

Remark 5.2. Lemma 5.1 extends to dimension n > 3 provided that for some e 
depending only on n and M 



sup 

Bi 



"("'^) -p{Q- 



< e 



sups^ |n(M,r)| 
for a three dimensional polynomial p G P2 and a rotation Q in M" . 

In the two-dimensional case the following lemma has been proved in [2, Corollary 
5.6]. 

Corollary 5.3. Let n = 3 and let u be a solution to (1.1) in Bi such that 
supg |u| < M and u(0) = |Vu(0)| = 0. Then there exist po > and rp > 
such that if 

sup |n(u,r)| > — 
Bi Po 
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for an r < ro then 



sup|n(u, 2 ^r)\> sup|n(u,r)| + JT]o/2 

Bi Bi 



TC \ \J ™^ '"' ^ ^^^^ ^osir/s) j forO < s < r, 



and 

snp\u\ > 

where 770 is the positive constant in Lemma 4-6. 

Furthermore, there exists a constant k ~ k(M, n) such that 

sup |n(u, 2~^r) I < sup |n(u, r) | + KJ, 

Bi Si 

sup|n(u,s)| < 2sup|n(M,r)| +Kj for se (2"(-''+iV, 2"-'V], 

Bi Bi 

and 

sup |u| < 16 1 ( - I sup \u\ + Ks^ log(r/s) ) . 
Bs \^rJ B^ J 

Proof. Lemma 5.1 applies, so that 

sup|n(u,2"V)| > sup|n(u,r)| +?/o/2 > — . 
Si Bi Pa 

It follows that Lemm.a 5.1 applies again with 2~^r. Thus wc may iterate and deduce 
that 

sup|n(u,2^^r)| > sup|n(u,r)| + j?/o/2, 

Bi Si 

which together with (5.1) proves the first part of the Corollary. 
We also notice that by Proposition 3.7 we have 

sup|n(u,r/2)| < sup|n(u(rx)/r2,l/2)-n(ii,r)|+sup|n(M,r)| < K + sup |n(u, r)|. 

B\ B\ B\ Bi 

Arguing as above we get 

sup |n(u, 2~^r) I < sup |n(u, r) | + jn, 

Bi Bi 

which together with (5.1) proves the second part of the Corollary. D 

6. Existence of a True Three-dimensional Singularity 

Corollary 6.1. There exists a solution u of (1.1) in Bi C M'^ such that 

lim ^(^^) = ^'+2/' _ ^2 
r-i-O sup^ |u| 2 

or 

lim "M ,_x^ + j^^ 

r-i-O SUPg^ \u\ 2 

Proof. The proof is similar to that of [3] , so we will only give a sketch. We define 
the operator T = T, : C"(B+) -^ C°'{B+) by 

AT{u)^-f,iu-u{0)) mB+, 

T{u) ^ M{^^^ - zA on dBi f^{z> 0}, and 

^=0 on{z = 0}nBi. 
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Moreover we impose that T{u) has cyhndrical symmetry, that is T{u){x,y, z) = 
gix^ + y^ , z) for some function g. The fmiction /e(i) is a smooth approximation of 
X{t>o} s-nd M is some large constant. 

By Schauder's fixed point theorem there exists an u^ such that T^(u^) ~ u^. We 
may pass to the hmit hm£_s.o u^ = u. Defining u(x) — u(x) — ""(0) for z > and 
w(x) = u{x,y, —z) — w(0) for z < 0, we see that u solves (1.1). From the boundary 
condition we infer as in [3] that ^^{r) < —M, which implies that sup^^^ |n(w, r)\ is 
also large. From Corollary 5.3 we conclude therefore that 

sup|u| >—((-) (sup|u| +77os^log(r/s) 



But then each limit of 



u(rx) 



( iiorigf'{bi'<:!»i.-o ' ^ ^<"" ""■ "-"^ ' <= '"•I)' 



sups^ \u\ 

as r — )■ must be a polynomial p e P2. Naturally, p will have the same cylindrical 
symmetry as u. Therefore p = {x^ + y'^)/2 — z"^ 01 p = z^ — (x^ + y^)/2. 
Last suppose towards a contradiction that there are two subsequences such that 
one converges to (x^ + j/^)/2 — z^ and the other to z^ — (a;^ +2/^)/2. By a continuity 
argument we obtain in this case a third subsequence and a limit that is neither 
(x^ + y^)/2 — z^ nor z^ — (cc^ + 2/^)/2, a contradiction. D 

7. Estimating w — n(ii) - Zn(„) 
The following Lemma is a direct consequence of Corollary 4.1 in [13]. 
Lemma 7.1. Let p be a second order polynomial in R" and ||p||l=°(Qi) = 1- T/ien 

^"({H < g}) 
iog(/:"({H < £}))!' 

/n particular, 

C^'HIpI < e}) < C{n, a)e" for every e G (0, 1) 
and all a < 1/2. 

The following Lemma is related to the two-dimensional result [2, Lemma 6.1]. 

Lemma 7.2. Let u solve (1.1) in Bi C M" such that sup^^ |u| < AL and u(0) 
|Vu(0)| = 0, and for some p < po and r < ro let 

sup |n(u, r)| > -. 
Bi P 

Furthermore let g^ he the solution of 

^9r = X{n(u,r)>0} ~ X{u(r-)>0} *^ -^1 , 

(7r = on dBi. 

Then for each a < 1/4, 

(i) ||Z?VlU^(Si) <^(M,n,a)|supBjn(u,r)||"". 
(ii) max(^supB^ |n(.gr, 1)|, sup^^ |n(5,., 1/2)|J < C(M,n, a) | sup^^ |n(u,r)|| 
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Proof. Let p = Il{u,r). Wc know that Agr = 1 when p > and u(rx) < 0, and 
that Agr = — 1 when p <Q and u(rx) > 0; in all other cases it is 0. By Proposition 
3.7 we also have that 

u{rx.) 



P 



<Cn 



Combining those properties we obtain that Ag,. ~ outside the set {\p\ < Cq}- 
From Lemma 7.1 it follows that 

WAgrWmB,) < (/:"({bl < C„W < C{M,n)\snp\p\r" for each a < 1/4. 

-Bi 

Standard L^-theory (see for example [21]) thus implies (i). 

Rotating and setting q := Il{gr,t) = J2^=i'^j^'ji where i = 1 or t = 1/2, we 
obtain 

P'gllL^CBi) < ^iII^VIIl^(Bi) < ^2| sup|n(u,r)||"" 

Bi 

and 

|aj|<C3|sup|n(7/,r)||"" 

Bi 

for every 1 < j < ri, proving (ii). D 

The following Lemina is related to the two-dimensional result [2, Lemma 4.3]. 

Corollary 7.3. Let u solve (1.1) in Bi C M" and assume that sup^^ \u\ < M, 
m(0) — |Vw(0)| — 0, and that for some p < po and r < tq, 

sup |n(M, r)| > -. 
Bi P 

Then 

sup \U{u, r/2) - U{u, r) - n(Zn(„,.), 1/2)| < C(M, n, a)(sup |n(u, r)|)-" 

Bi Bi 

for each a < 1/4. 
Proof. For each r write 

u{r:x.)/r^ = n(u, r) + Zu{u,r) + 9r + K 

where Acjr = A{u{rx)/r^) - AZn(„,r), A/i,. = and gr{0) = |V.gr(0)| = /ir(0) = 
|V/ir(0)| = \U{gr, 1)1 = |n(/i,., 1)1 = 0. Next denote by gr the solution to 

Agr = Ag,. in Bi , 
(7r = on dBi. 

Then gr + hr = gr + hr for some harmonic function hr in Bi. From Lemma 7.2 it 
follows that 

sup|n(5,,l/2)|<Ci(M,n,a)(sup|n(«,r)r", 

Bi Bi 

and from Lemma 3.6 we infer that 

sup \U{hr, 1/2)1 = sup \n{hr, 1)1 - sup |n(.g,, 1) - n{gr, 1)1 

_Bi Si Bi 

< Ci(M,n,a)(sup|n(u,r)|)-". 

Bi 
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Thus in Bi 

Uiu, r/2) = niu, r) + n(Zn(„,,), 1/2) + Uig,, 1/2) + n(/i„ 1/2) 

and 

|n(u, r/2) - n(w, r) - n(Zn(„..), 1/2)1 < C(M, n, a)(sup {Uiu, r)|)"". D 

8. Classification of Blow-up Limits in M^ 
— An Unexpected Symmetrization Effect 

In this section we will show that if lim.-^oo " ^^^i i = v, where b is a harmonic 

rj 

polynomial, then p = 2xz, p = (x^ + y^)/2 — z^ or p = z^ — (x^ + y^)/2 up to a 
rotation. 

Theorem 8.1. Let n = S, let Au = — X{u>o} *'^ Bi and assume that u{0) = 
|Vu(0)| = and that the monotonicity energy satisfies lirur-i-o 'i'x"('') ~ ~oo- Then 
each limit of 

^^-n(u(x),r), 

as r -^ 0, is contained in 

{Z(gx) : g e 2} U {^i(gx) : g e 2} U {^2(gx) : g e 2}; 

/lere 2 is the set of all rotations of M.^ . 

Proof. Suppose towards a contradiction that the statement is not true. By Propo- 
sition 3.2 (i) there exists a solution u and a sequence rj — ^ such that after rotation 



.1) lim 

j->-oo 



w(rjx) 



suPr . \u\ 



-PSo\ =0 



for some So e (0, 1/2) and ps„ = (1/2 + So)x^ + (1/2 - 6o)y^ - z^ and ^o e (0, 1/2) 
or psa ~ z^ — (1/2 + So)x^ — (1/2 — <5o)2/^. We may assume that ps„ = (1/2 + 
5o)x^ + (1/2 — So)y'^ — z^ . Furthermore, from Proposition 3.2 (i) and Proposition 
3.7, limj._j.o sup^^ |u(rx)/r^| = oo. 

We are going to prove a decay estimate for 5^ = 5^{u{r-)) in r which will lead 
to a contradiction to (8.1). 

By Theorem 4.5 (vii), 

3A„((5) - A{5) 1-25 ^ , , 
^^^^l^:^!^ >^^ for. e (0,1/2). 

Thus 

where a; is a continuous function on [0, 1/2]. 

By Corollary 7.3, using Corollary 5.3 to estimate 

(sup|n(u,r)r"=0(|log(r)|-"), 

Bi 

we obtain for every a < 1/4 that in Bi, up to a rotation depending on r, 

(8.3) |n(u, r/2) - TrPsA - n(ZxA, 1/2)1 < 0(| log(r)r"); 
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from here on, Z^ is the unique solution to 

AZa = -X{pi>o} in M^ 

Z5(0) = |VZ5(0)| = n(Z,, 1) = lim,,,!^^ ^5(x)/|x|3 =. 0. 

In particular, for the Kq defined in (4.1), 

n(Z,, 1/2) = -Ko[(3A,(J) - A(J))a;2 + (zAy{S) - A(<5))y2 + (3A,(5) - A(<5))z2 

= -KoA'i[(l + 18)x^ + /r2(l + 2(5)y2 + ^3(1 + 2<5)z2], 

where — using the fact that 11(^5, 1/2) is harmonic — 





(l + 2(5o) ' 






3Ay(<5) - A{6) 
' 3yl,((5) - A(,5) 


1-25 + k((5) 
1 + 2(5 


and 








3A.(<5)-A(<5) 


-(1+^2) = 


It follows that 







^(5) 



1 + 2(5 



n(Z5-, 1/2) = -Ko/-i:i[(l + 2(5)x2 + (1 - 2(5)y2 + ^((5)^2 - 2z^ - ^((5)2^] 
= C5(2p5 + K((5)(2;2-z2)) 

for a = -Ko(3A:,((5) - A((5))/(l + 2(5) > <: > 0, 5 e [0, 1/2] (see Theorem 4.5). 
Invoking (8.3), this implies that in B\. up to a rotation depending on r, 

(8.4) |n(u, r/2) - r,v,A - c,a (2p,.4 + K[8^){y' - z")) \ = 0{\ log(r)|-"). 

The fact that c^ > c > as well as the estimate k > a; > consequently prove 
together with Corollary 5.3 that, rotating slightly. 






(8.5) (5,^/2 < K - cinrfr + C2I iog(r)| 



Note that estimate (8.5) is independent of rotations. As long as a;((5^_fc^ ) > 
I log(2~''ro)|^"/^, induction of estimate (8.5) in k yields a logarithmic decay of 8^ 
in r. It follows that J^ -^ as r — > 0, contradicting the assumption limj_>.oo ^r ^ 
(5o > 0. ' D 

9. Unique Tangent Cones at True Three Dimensional Singularities 

From the previous section we may infer by a continuity argument that in three 
dimensions, assuming w(0) = |Vm| = as well as lim,,_j.o ^"^air) ~ — oo, then one of 
the following three statements holds: 

(i) linr M'Q^p^) _ n(., r)(Q(r)x)) = Z(x), 

(ii) / »(rQ(r)x) _n(^,,)(Q(,)x)) =Zi(x), 

(iii) Una /'!^(!:Q(!M _ n(^., r)(Q(r)x)) = Z2(x) 
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for some Q{r) G Q. However at this point we do not yet know whether the rotation 
Q{r) converges as r — > 0. 

In this section we will show that in the case (ii) and (iii), "^7^-^ — n(u,r)(x) 
converges as r — > 0. In Section 11 we will show a similar result in the case (i). 

Theorem 9.1. Let n = 3 and let u solve Au = ^X{u>o} *'^ ^i such that m(0) = 
|Vu(0)| = and M := sup^^ |u| < +00. There exist constants r{M) > 0,c(M) > 
and K{M) < +00 such that if 

s e (0, r(A/)), n(w, s) > K{M) and 5^(u(.sx)) < c{M), 

then there is a rotation Q such that either 

u{rx) 



(9.1) 



-n(u,r)(x) ^ Zi(Qx) asr^O 



(9.2) 



u(rx) 



2 - n(u, r)(x) ^ Z2(Qx) as r ^ 0. 
Moreover, there exist /3 > and C{M, /3) < +00 smc/i t/iai 
n("'^) ^ <C(Af,/3)(if(Af) 



sup 



suPsjn(M,r)| supbJpI 

/or all r g (0,s); /lere p{x) = {x^ + y^)/2 
-(a;^ + j/2)/2 - z^ in the case (9.2). 



loe 



z^ in the case (9.1) and p(x.) 



Proof. Observe that the assumptions imply by Corollary 5.3 as in the proof of 
Corollary 6.1 that r^ > K{M) for r < s and that r^ — >■ +00 as r ^^ 0. Moreover 



we see from Theorem 8.1 that 



'u(sx) 



— n(u, s)(x) is after rotation close to Zi(Qx) 



or Z2(Qx). We may assume that it is close to Zi{x). 

We will follow the strategy explained in the proof of Theorem 8.1, and use the 
notation of that proof. Remember that by (8.4) and Corollary 5.3, up to a rotation 
depending on r, 

rPsA - csa{2psa + n{5^){v^ - z^))\ = 0{t-''), 



(9.3) |n(u,r/2)- 

where c^ > c > and 



(9.4) 



k{5) = (1 + 2J) 



iAy{5) - A{5) 



(1 - 25) > uj{S) > 0. 



'3A,i5)^A{5) 

In Theorem 8.1 we worked to exclude the case that 6^{u{r-)) S [f3,l/2 — /3] for 
positive /3 and small r, and in that (5-rcgime, uj has been bounded from below by 
a positive constant. In the present proof, however, we are interested in the regime 
(5 — >■ 0, where uj degenerates. In order to deal with this difficulty, we will make a 
Taylor expansion of 3Ax{5) — A{6) and 3 Ay {6) — A{S) at the point 6 = 0: From 
Theorem 4.5 wc infer that 

(9.5) 3A^{6) - A{S) = 3A^{0) - A{0) ^S + 0{6^) 



and 
(9.6) 



3V3 



Att 



3Ay{S) - A{S) = 3Ay{0) - A{0) + -^S + 0{S') 

3v3 
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Plugging this information into (9.4), we obtain that 

(9.7) 4S<K{S)<(^5+^y. 

Dividing (9.3) by r^, rotating slightly and recalling that ps = (1/2 + 5)x'^ + (1/2 
S)y'^ — z^ and using that q > c > we infer that 



(9.8) ^t^<^^~c^— + C2T~^-'^ 



where c\ > is a universal constant and C2 < +00 depends only on M and a. 
Note that estimate (9.8) is independent of rotations. In the following three Claims 
we will describe how (9.8) leads to a decay estimate for 8^ . 

Claim 1: There is a universal constant (3 > and C3 = C3{AI) < +00 such that 
ifTr > C{M) and 5^ < r"^, then 

Proof of Claim 1: First, (9.8) as well as the assumption in the Claim imply that 

S^/, < (1 - (ci - C2r-"-^)r-i)r-^. 
On the other hand, < Tr/2 — Tr < k (see Corollary 5.3), so that 

(9.9) (^lL^y^>i-p^r-\ 

It follows that, provided that /3 has been chosen small enough (depending only on 
the universal constants k and ci) and r^ is large enough (depending on k,ci and 
M), then 



s^/2 < r;f,, 



proving the claim. 



Next we consider the case when 5:^ > t,. ^ . 
Claim 2: There is a constant /3 = /3(a) G (0,a) sucli tliat if 

<5;^<C4(M), Tr>C5andS^>T-^, 

then 

(9.10) s^^,<{l-^)sf. 

Moreover, ifS^^^,. > '''2-''r ^'-"^ each k < ka then 

,A ^rfS^ 1 



2^ ''Or 2^'=0r 

Proof of Claim 2: Equation (9.10) is a direct consequence of equation (9.8). The 
last part of the Claim follows from an induction of the first part (noting that the 
assumption 5^ < C4(M) is satisfied inductively): If S^-k^ > '^2 ^r ^'^^ each k < ko 
then 






Cl 
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That product can be estimated for Tr > C5, calculating 

fcfj — 1 fco— 1 



iog(ntro^( 



1 



2t,- 



2-fcr 



fc=0 
fco-1 

fe=0 

where we have used Corollary 5.3. Thus 

Cl 



2T2-kj. 



Cl 



k=0 



4rn-fc 



2-fcr 



Cl < '^1 1 . / fepK + T,. 

4(Tr + kA:) ~ 4 °^ V Tr 






< 



=1/4 



2T2-kj,/ VfcgK + Tr 

Choosing /3 even smaller such that 2/3 < ci/4 and using once more Corollary 5.3, 
we obtain 



"2 "«'■- VfcoK + r^ 

1 
< 



/^ T^^d^ 



< 



T^^S^ 1 



^ 2^''0r 2^'=0r 



proving the Claim. 



Claim 3: There is a constant (3 = P{a) G (0, a) such that if T2-fcor ^ C(Af) and 
S^ < Ci{M), then for k>ko, 

<5^-fc„ < T~_^ for k > ki and 



J2-fcr — '2-^'r 
02-fcr — 2/3 



for k < ki 



'2-i'r 

for some fci < — T^{d:^)'e . 
Vo 

Proof of Claim 3: We apply Claim 2 up to the first fci such that S:^^^^ < '''o-''! ' 
and we apply Claim 1 for k > ki. From Claim 2 and Corollary 5.3 we infer that 

fci<-r2(5,^)i. 

Observing that the assumptions for T2-k^ are satisfied for k > ko hy Corollary 5.3 
finishes the proof of Claim 3. 



In the last part of our proof we will use the decay estimate in Claim 3 in order 
to estimate how much n(u(r-)) moves when varying r. Let r^ := 2~ s, r^ := T2-kg 
and 6^ :~ 5^-k • First, we infer from (9.3) that up to a rotation depending on fc. 



(9.11) sup 

-Bi 



n(u,rfe) \l{u,rk+i] 



< sup 



Tfc Tfc+i 

TkPs^ TkPsA + CgA {2pgA + K{5^){y^ - z^)) 



Tfe 



Tk+1 



+ Ci{M,a)T^ 



(l+a) 



where < c,. = -Ko(3A,(5,^) - A{5t))l{l + 25^) < Cg, 4<5^ < niS^) < Cj6^ 
and Cg, C7 are universal constants. Another fact we infer from (9.3) is that 

(9.12) Tk+i=Tk + 2csA+0{5^). 
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Plugging (9.12) into (9.11) yields 

n(u,rfe) Il{u,rh+i) 



sup 



< sup 

Bi 



Tk Tfc+i 

-Aa. + TkOisi) 



+ Ci{M, a)r-(i+"^ < C8(M, a)T-'^+") + Q 



rfe(rfc + 2c,.+0(<5,<l)) 
where Cg is a universal constant. Iterating this estimate we obtain 

n(u,rfe) \l{u,rk+m) 



Si 

Tk 



(9.13) 



sup 

Bi 



Tfc 



Tk+r 






From Claim 3 (applied twice) and Corollary 5.3 we conclude that, choosing riM) 
small enough such that t^(j\/) > CiJVl')^ setting -p = x^ ^ y^ — 2z^ and letting 



■nij - 


-> oo, 






sup 

Bi 


n(u,rfc) 

Tfc 



< Cio(M, /3) 5: r-^-'' + CrfSt Y. r-'' 



2/3 



i^k 



< CMM, P) Y,{M + ^77o/2)-l-'^ + Cisrf <S,^ ^(M + zr;o/2)-i 



-1-2/5 



i—k 



i—k 



< Ci3(A/,/3)t7^ for all fc > k{M). 
Using once more Corollary 5.3 we obtain the estimate of the Theorem as well as 

lim /^ ^(^"+^^) _ n(M,r,x°)(x)^ = Zi(x) for x" = 0. 



D 



Corollary 9.2. Lef n = 3, /ei u solve (1.1) in Bi and suppose that 

(9.14) hm ^^i!ZE) ^ (2:2 ^ 2)/2 _ ^2 ^^ ^j^^ ^i(rx) ^ ^2 _ ^^2 ^ 2^/2. 

r-^O SUps^ \u\ r^O SUpg^ |u| 

T/ien there exists an rg — ro(u) and f,g G C'^'^(i?' ) such that 

Br„r\{u = 0} = ({(a;,y,/(x,y)) : {x,y) e B;Ju{(x, y,.g(x, y)) : (x,y) e S;j)nB, 



Moreover /(.T, y) - Va?MV/V2 e C1(B; J and5(x,2/) + V^^^lV/%/2 G C1(B;J. 
The Lipschitz- and C^ -norms corresponding to the above statements are uniformly 
bounded for solutions v sufficiently close to the fixed solution u in L°°{Bi), provided 
that each v satisfies 

u(rx) 



hm — - 



lim 



r-i-O SUp^^/^u") \v\ r-s.0 SUp^^ \u\ 

for some rotation Q'" at a singular point ^^ sufficiently close to the origin. 

Proof. We will show that {u = 0} n S+ = {(x, y, f[x, y)) : [x, y) € S;j n Br,, for 
some / £ C°'^(i3^|J and /— \JliP^^\^]p- j \f2 G C^. By symmetry a similar statement 
holds in B~ . We will also assume, for the sake of definiteness, that 



Wr(x) 



u(rx) 



^ [x^ + y^)/2 - z^ in C^-''{Bi). 
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By the C^ '"-convergence, 



and 



sup 

Bi 


It-- 


< ^ 


sup 

Bi 




2 



<uj{r) 



for some modulus of continuity aj(r) — > as r -> 0. It follows that {v^ = 0} n {Bi \ 



B1/2) C {(x, y, z)£Bi\ Bi/2 ■■ dist(-, {a 
of continuity a. Therefore 



2-z }) l£ "'('')} for some modulus 



dVr 1 -I , „ 

-^ <-- on B+\B,^,. 

From the implicit function theorem and C^'"-regularity we infer that {vr = 0} H 
{B^\Bi/2) is a C^'"-graph with bounded C^'"-norm (independent of r). It follows 
that {u = 0} is the graph of a Lipschitz function / in B^ and we only need to 
show that f{x, y) - yf^^TWl^ € C^{Br„). 

We know that / £ C^'°'{B'^\^B'^ for every s > and that /(rx, ry)/r is bounded 
in C^'"(i3j^ \ ^5./2)- Thus it is sufficient to show that 



lim |V(%/2/(a;,y)- v/a;2+y2)|^o. 

(x,y)-)-(0,0) 



Let us consider any sequence {xj ,yj) — > and denote 



Tj . Then 



f{rjX,rjy)/rj will converge to \/x'^ + y'^/V^ in C^'^{B[ \ S1/2), implying that 
(v/(a:,2/) 



.\/xFTy^ 



V2 



{x,y) = (xj,yj) 



As the sequence {xj,yj) is arbitrary, it follows that / — i/a;^ + y'^/^/2 <E C^ . The 
uniformity follows from the uniformity in Theorem 9.1. D 



10. Stable Cones 
Theorem 10.1. In M"^ there exists a solution of (1.1) in Bi such that 



(10.1) 



0<l f \Vw\'< f \Vw\-' 

^ JBi JBi 

^1,2 



Bin{u=o} 



w 
Wu\ 



dU' 



{or each w e Wg' (Si). Moreover u (f. C^'^(Bi/2). 

Notice that the right-hand side in (10.1) is the second variation of the energy 
/b (I^'^P/^ - max(it,0)) of equation (1.1). 

Proof. By Corollary 6.1 there exists a solution v of Au = -~X{v>o} in -Si such that 
the blow-up limit at the origin is Zi. 
Let 

^(sx) 
"W := — ir- 



for some small but fixed s. 
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For some large M to be determined later and sufficiently small s, Theorem 9.1 
together with Corollary 5.3 yields that 

(10.2) |Vu(x)|>(Af + log(^))|x| 

on r = Bi n {u = 0}. 

Choosing s if necessary even smaller, Corollary 9.2 implies that F consists of two 
Lipschitz graphs in Bi. 

Note that since the origin has zero capacity we may by a limiting argument 
deduce that the second variation is well defined for aU w € Wq' {Bi). 

If w £ W^'^ then w\r & VF^/^'^(r) by the trace theorem, which is valid for our 
Lipschitz free boundary. Also from the trace theorem, combined with Poincare's 
inequality, wc infer that for each w € Wq' (Bi) 

\\w\\ivi/2.2^-P) < C'i||w||h'i>2(Bi) < C'2||Vw||l2(Bj). 

Using the Sobolev embedding, we obtain for w G Wq'^{Bi) 

(10.3) l|w||L-i(r) < C'3lkl!vKi/2,2(r) < C4|lVu;||i2(s,). 
Thus, using (10.2) and (10.3), 

\ 1/2 / /. \ 1/2 

dn' < { I 7^-— T dn'^] / w'^ dn^ 



IVul - V ./r IVmP 



r |VU| \Jy I v'U|- / \Jy 

1 \ 1/2 

r(A/+(log(rx|-i))2)|xP / JB, 
On the other hand, |F n dBr\ < C^r by Corollary 9.2, so that 



<cA I .,, ,. ,:, ,,,„., ,„ dHM / \vw? 



^n^"',cjr -,,... )-^^ 



,r(A/ + (log(|x|-i))2)|x|2 ; - 'V7o (M+(log|r|)2)r2; " ^m' 

Choosing M > AC^C^ we arrive at 



r |Vu| - 2 jB^ 



D 



11. Unique Tangent Cones at Unstable Codimension 2 Singularities 
In Theorem 9.1 we showed that if 

lim ^^"^-^^ ={x^ + y^)l2-z\ 
J->oo sup^ |m| 

^ j 

then the blow-up limit is unique and wc obtain a quantitative convergence estimate. 
In this section wc will show the corresponding result in the case that 

lim \ , = x^ - z^. 

j->-oo supg \u\ 

This case corresponds to (5 = 1/2 in the notation of the previous sections. To 
make Taylor expansions of 'iAx{5) — A{5) etc. around the point 5 = 1/2 would 
be rather clumsy. To get around that we will change the parametrization to ps = 
(1 — 5)x'^ + 5y'^ — z^ and use the B^, By, B^^B defined in Section 4. 
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Theorem 11.1. Let n = 3, let u solve Au = — X{u>o} ™ Bi and suppose that 
M :~ swpB W\ < +oo,x'^ G _Bx/2,w(x'^) = |Vu(x*')| = and that there exists a 
sequence Vj — > such that 



lim 



where Q-^o is a rotation depending on the point x*' . Then the limit 
lim ^ ^ -n(u,r,xO)(x) = Z(QxOx) 



exists (and is thus unique). 

Moreover, for each 7 € (0,1/4) there exist constants r{M,^) > 0,c{M) > 
0,K{M) < +00 and C{M,j) < +00 such that 

s e (0,r(A/,7)), U{u,s) > K{M) and 5^(w(sx)) < c{M) 

imply that 

n(u,r,x")(x) 



2(Qx«x)i(QxOx)s 



< 



C(Af,7)(AlM) + |log(^)|)" 



supsi |n(u,r,xO) 
for all r G (0, s). 
Proof. We may assume that x" = and that the rotation Qx" is such that 

2(QxOx)i(gxOx)3 = Xi - X3. 

By CoroUary 7.3, up to a rotation depending on ?■, 
(11.1) sup |n(w, r/2) - T.p,., -n(Z 1/2)1 < C(A/,n,7)T-T 

for 7 < 1/4 and r < r{M). FoUowing the strategy in the proof of Theorem 9.1, we 
arc going to use (11.1) together with an analysis of n(Zp ^ , 1/2) to derive a decay 

estimate for (5^(u(rx)) > (cf. Definition 4.2) in r. That decay in turn will make 
it possible to estimate how much n(w(r-)) moves when decreasing r. Note however 
that as the singularity examined in the present section is by [15] unstable, we cannot 
expect to obtain the decay by a simple iteration as in the proof of Theorem 9.1. 
The "pinning effect" of the convergence assumption 



lim 



^ ^2 -n(^,r„x°)(x)j =Z(QxOx), 



has to enter the proof. 
Claim: For r < s, 

2C(M,7)t-'' 



(5f < 



KoiSCyid^) ~ Coi6^)) 



Proof of the Claim: As the proof will be concluded by a continuity argument in r, 
we assume that U{u,r) > K{M) and 5^{u{s:x.)) < c{M). From (11.1), (4.1) and 
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Theorem 4.7 we infer that in Bi, up to a rotation, 



Ii{u, r/2) = TrPsB + n(Zp „ , 1/2) + 0{t--') 



(11.2) = [Tril - 5?) + Ko((3B,(l/2 -5)- 5(1/2 - 5))x'' + (3B,,(l/2 - 5) 

- 5(1/2 - 5))y^ + (35,(1/2 - 5) - 5(1/2 - (J))^^) + Oir--^) 

= {Tril - 5f ) + Ko(l - Co(5f )d:f + o{5f)))x^ 
+ (r,5f + Ko(3C,(,5f )<5f - Co(5f )(5f ))2/2 
+ ( - T, + Ko(-3Cy((5,^)(5f + 2Co((5f )5,f - 1 - o{5f)))z'' + 0(r-^). 

Rotating the coordinate system shghtly to x,y,z, we deduce from (11.2) that 
the quotient of the y^ and the a;^ coefficient of Il{u, r/2) is estimated from below 

by 

^f/2 ^ T,Jf + 3KoCy{S^)S^ - KoCn{S^)S^ - Cit^' 



- V2 Tr{l-SB) + Ko-KoiCo{6^)d^+o{SB)) + CiT,r^ 
We maintain that for S^ < C2(Af, 7) and t^ > C3(M, 7), 

Subtracting the two quotients we end up with 

■ l-^f/2 1-^.^ 

> s[Ko{3Cy{6^) - Co{6^))S^ 

- KoiSCyiS^) - 2Coi5^)){S^f - Ko(5f + 5f o(,5,f ) - Ct-^] , 

where 

S = ((r,(l - 5^) + Ko - Ko(Co(<5f )5f + o(<5f )) + Ct-^)(1 - 5f ))"' . 

For D to be non-negative — 6^ being by assumption small, r^ being large and 
Cy{S^) sa Co((5f ) « I log(5f I by Theorem 4.7— it is sufficient that 

(11 4) J^ > '^^IZ^ 

Thus (11.3) holds, and 

^'^^^'' ^ < <5f < C2(Af, 7) implies that ^^ > ^'' 



Koi3CyiSB) - CoiSB)) - "^ - ' ' l-Jf/2 - 1-5,^' 

which in turn implies that S^,^ > 5^ . Using Theorem 4.5 (i) together with Theorem 
4.7 (i) while observingthat -§^{5^) = -j^{k~^?) and that |^(<5f ) = -§^(5- 
(5^)), we conclude that 

{3Cy{5?,2) - Co(<5f/2))'5f/2 > (3C,(5f ) - C,{5f))5f. 
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As T^l < Tj. ^ , it follows then that 



(11.5) 



<55.. > 



2^1-72 



'■/2- Ko{3Cy{S^)-Co{5n)' 



In this case S^^^^ ^ ^r for fc = 1, 2, . . . . Altogether we obtain that 



either limini S^-k^ > oi S^ < 



2CiT-^ 



Ko{iCy{SB)-Co{S^))' 



Since lim inf fe_j.oo <^2-'=r ^ ^ would contradict our assumption that 



n(u, Vj)/ sup |n(M, rj)\ -^ x^ — z^ as j — > oo, 



we have proved the Claim. 

In the last part of our proof we will use the decay estimate in the Claim in order 
to estimate how much n(u(r-)) moves when varying r in the interval (0, s). First 
note that the Claim and the fact that Cy w Co when i5 << 1 imply that 

(11.6) Cy{5?)5f < C4T-\ 
Next observe that by (11.2) and (11.6), 

(11.7) T,/2 =Tr+Ko+ OiCyi6^)6^) =Tr+Ko+ O(t-T). 

Using (11.2) once more along with (11.7) and (11-6) we obtain 
Il{u,r) n(M,r/2) 



sup 

Bi 



< sup 

Bi 



supBjn(u,r)| supBjn(u,r/2)| 



JTrjl - S^) + Kq - Ko{Co{5^)S^ + o{5^j) + 0(t-^))x^ 
{Tr6? + 3KoCy{S^)S^ - KoCo{6^)S^ + 0(T-^))y^ 

Tr+Ko + 0{T,r'^) 

( - T, + Koi-3Cy{6^)S^ + 2Co(<5f )<5f - 1 - o{S^)))z^ 



Tr + Ko + 0{Tr ^) 

As in the proof of Theorem 9.1, an iteration leads to 

n(M, 2-'=s) 



<C5rr(i+^). 



supsi |n(M, 2 ''s 
and we obtain the desired estimate as well as 

'ii(x'^ + ?'X 



{x' - z') 



< Cqt Jt, 



lim 

r-s-O 



n(M,r,x")(x) = Z{Q^o^). 



D 
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12. Structure of the Singular Set in M^ 

So far wc have shown that if Au = —X{u>o} in Bi C R'^ then the singular set 
S'" = {a; G -Bi : m(x) = |Vw(x)| ~ and hmr-yo 'i'x(^) ~ — oo} is divided into two 
parts Sf = {x e 5" : hm.^o "^[.g^'p = ^ - z" for some Q G Q} U {x e 

5" : hm.^o tup^X? = -(^ - ^') fo^' ^°'™ Q ^ 2} and ^2" = {x e 5" : 
hnir_j.o su ImI ^ ^"^ ^°^ some Q £ Q}. In this section we show that Si consists 
only of isolated points, and that S2 is locally contained in a C^-curve. We also 
derive a compactness result for 5*2 • 

Lemma 12.1. Let n = 3, let u solve (1.1) and let x*^ e 5". Then there exists an 
r = r(w,x") > such that {u = 0} n {Vw = 0} n Br(x") = {x"}, that is x" is the 
only singular point in a small neighbourhood of x*^ . For each class of solutions v 
sufficiently close to u in L°°{Bi), {v = 0} Ci {Vw = 0} contains at most one point 
in Br- 

Proof. Suppose towards a contradiction that there exists a sequence of solutions 
u^ — > u in L°°{Bi) as well as sequences {u^ = 0} n {Vu-' = 0} 9 x^ -^ x" and 
{u3 = 0} n {Vu-?' = 0} \ {x-J} 9 y-J' -)■ x° as j -^ 00. Let r^ = |yJ - x^]. Then, 
passing if necessary to a subsequence, (y^ — 'xP)/rj ^- ^ G dBi. On the other 
hand, by M^^.p.j-egularity of the solution, u^ -> u in L°°(Bi) n W'^^p{Bi/2) so that 
the assumptions in Theorem 9.1 are satisfied in Bp(xP) for small p and sufficiently 
large j. Rotating each solution suitably around the origin, we obtain that x-' G S^^ 
and that 

(12.1) lim ^' I ,/ =w{^)^±{^^-z'). 

j->oo sup^ Im-^I 2 

By C^'"-convergence in equation (12.1) it follows that w{^) = |Vw(4)| — which 
is a contradiction since |^| = 1 and the origin is the only point where w — \Vw\ = 
0. D 

We continue this section with a regularity result for 6*2 • 

Theorem 12.2. Let n — 3. // G 5*2 then there exists an r{u) > such that 
S'2 n B^ty^\ is contained in a C'^ -curve. For each class of solutions v sufficiently 
close to u in L°°{Bi), the curves containing 5*1 are relatively compact in C^{Br(u))- 

Proof. Let us consider a sequence of solutions w-' — > u in L°°{Bi). By uniform 
M^^'P-regularity of the solution, for sufficiently small s > 

n(u,s) > 2K{M) and S^{u{sx)) < c(Af)/2, 

and for all sufficiently small |x"| and all sufficiently large j, 

n(uJ,s,x°) > K{M) and (5^(u(x° + sx)) < c{M). 

From Theorem 11.1 we obtain therefore that 



sup 



u^ (y + rx) ■ 
TT -Pi 



<ei 



supsi |uJ(y + ?■■ 

for all sufficiently large j, all y G ^3 n Bp, a rotation Q^, p{xi,X2,X3) = 2x12:3, 

p^(x) =p{Q^yX) and ah r G (0,ri). 

Uniform cone flatness: For each e > there exists an Sf > such that for 
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sufEciently large j and all y e Sf n Bp,, {u^ = 0} n {\\7u^\ = 0} n Bs{y) C 
{y + sQy{xi,X2,X3) : xj+xl < exl} for s e (0, Sc). 

Proof of uniform cone flatness: Suppose towards a contradiction that there exists 
an eo > 0, a subsequence of solutions 



—^ ' y^ ' ^ 2x1X3 in C''"iB,) 



sups^ liPiy^ +Sj-)\ 



and a sequence of points $,'' ^- $, £ dBi such that v^ {$,■') = |Vw-'(^-')| ~ and 
(^j )^ + (^i? > eo. Then (^0)2 + (^0)2 > eo, contradicting = |Vi;0(|")| = 
2-\/(fi)^ + (^3)^ ^i^d thereby proving uniform cone flatness. 

A standard consequence of the uniform cone flatness is that the class of curves 
containing S2 H Bp^ is for large j relatively compact in C^. An argument by 
contradiction yields the Theorem. D 

The following corollary can be regarded as an extension of [15, Corollary 7.2] outside 
a small cone (even outside a cusp) in the y-direction. 

Corollary 12.3. Let n = 3 and suppose that for some solution u to equation (1-1), 



u{rjx) 

lim , , = 2xz. 

j->c» sup^ \u\ 



Then for each > there exists an r{u) > such that {u = 0} n Br[u) H {|yp < 
0(|x|2 + I^P)} consists of two 2-dimensional C^ -manifolds restricted to B^^u) H 
{|yp < 0(|xp + |zp)}, intersecting at right angles at the origin in the xz -plane. 
For each class of solutions v sufficiently close to u in L°°(Bi) and having each an 
S2 -point sufficiently close to 0, the manifolds are relatively compact in C^ . 

Proof. The proof is similar to the proof of Corollary 9.2 and left to the reader. D 



13. Appendix 

Proof of Theorem 4-5: 

Let us begin by proving (vi), that is 

(13.1) 

3(4'(5)-<(5)) + 2<S(34'(<5) + 3A^(<5)-2.4"(<5)) + 2(34(<5) + 3A;(<5)-2A'(<5))>0 



for S e (0, 1/2) 
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For (5 > we have that jVp^l 7^ on {ps = OjnQBi and we may thus differentiate 

, Ay 



Ax, Ay and A. Thus 



r/^ ., , ^- J -, — r,, n, ,,aarccot(v'l/2 + (5cos(2(/))) 

= 8 / sin^ (arccot ( ^1/2 + ^cos(20)) ) cos^ ((/>) ^^ V^ !-— ^ < 

= 8 /" ^ (^ ( cos(3arccot(Vl/2 + <5cos(20))) 



- 9 cos(arccot( v/l/2 + 8 cos(2(?!))))) j cos^(0)d0 

32 r/2 5 .(2 + ,5cos(20))v/1 + 2(5cos(20)n 2/.^,. 
= -y X 5^ V (3 + 2^cos(20))3/2 ) ''' (<^)^<^' 

= 8 /" — (^ ( cos(3arccot( v/1/2 + <5 cos(2(/.))) 

-9cos(arccot(v/l/2 + (5cos(20))))') sin2((/))d(/) 

32 r/2 g /(2 + <Scos(20))v/l + 2(5cos(20)x . 2,,^,, 
^^Yio a^V (3 + 2^cos(20))3/^ jsm(^)d^ 

and 

A' (,5) = -8 / — (cos(arccot(\/l/2 + ^cos(2</)))))d0. 
Jo c/" 

Differentiating once more, 

,„,,,_ 32 /-"/^ 92 ^ (2 + ^cos(20))v/l + 2^cos(2(/.) ^ 2 
^^'^" 3 /n 9<52V f3 + 2Jcos(2<i<))3/2 ;c°^ 



32 r/2 92 ^ (2 + ^cos(20))Vl + 2Jcos(20) ^ . , 



and 



3 Jo dS^\ (3 + 2Jcos(20))3/2 

32 r/2 92 /(2 + ^cos(20))^l + 2,5cQ 

yjo 5^^ (3 + 2<5cos(20))3/2 

^/2 ^2 



A"(5) = -8 / —2 ( cos(arccot(^l/2 + 5 cos(2(f>)))) d(f>. 

Using the last three identities we may write the left hand side in equation (13.1) 
as the sum of the following three terms (13.2), (13.3) and (13.4): 



''/^ 52 ^ (2 + ^ cos(20))^l + 2<5cos(20) 
(3 + 25cos(29i))3/2 
/2 92 ^l + 2<5cos(20) 



(13.2) 3{A;-A'^)^32J^ — ( 

(13.3) 2.(3A^' + 3A; - 2A") = -32S f" ^ ( ,f^^^7,tt ) 
^ ^ ^ ^ ' Jq 9(^2 V (3 + 2(5 cos(20))3/2y 



cos(2t/))(i0. 



and 
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In order to estimate (13.2), (13.3) and (13.4), we will the change of variables x = 
(5cos(2(/)) and then use Taylor expansions: First we notice that — using double 
factorials {2k + 1)!! = 3 x 5 x • • • x (2A: - 1) x (2/c + 1) — 



yrT2^ = 1 + :, + g(-l)^-+i (^^Jll:, 



fe=2 

and 



fc 

fc! 



1 - ^ V^ .^.(2fc+l)^. 

fc=0 



(3 + 2.t)3/2 3^/3^' ' 3^fc! 
both sums are absolutely convergent for |a;| < 1/2. Thus 



X , 



VT+2x l + xf^^ _,,j2fc + l)!! 



(3 + 2x)3/2 3V3 



E(-)'^^I4F^') 



fc=0 



J2fc+1)!!^, 



fc=2 ^ fc=0 

1 , 1 V^. ,,,+i/(2fc-l)!! 



3V3 3V3^' ' V 3*^-fc! 



^D-i)'«(^^l^(->))^ 

J2fc+1)!! , 



A;! ^^3V3fr^' ' 3^fc! 

Notice that the product of the last two sums will equal a sum X^fcLo '^kx'', where 
Ofc > when k is even and a^ < when k is odd. Inserting a; = ^ cos(2(/)) and these 
Taylor expansions in equation (13.3) and using that a^ > for even k and that 

Jq Ofc cos'^ (2 (f>)d(f> ~ for odd fc, we see that 



2(5(3A;' + 3A;' - 2A") 



Similarly we may estimate the left-hand side in (13.4) and obtain that 

2(3A; + 3^;-2A')>0. 

Next, we make a Taylor expansion of the integrand in equation (13.2). First, we 
calculate 

{2 + x)VlT2i^{2 + x)(l + x + Y,{-lf+'^^^^^^p^x'') 

k=2 

= 2 + 3. + £(-1)^+^ (^^^Tr^(3fc - &))xK 

k=3 
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Therefore 

(2 + a;)Vl + 2a; 



(3 + 2x)3/2 



3V3 



fc=0 



fe=3 ^ fc=0 



3^/3V ^1^ ' V 3fefc! 



2 + X + ^(-1)^+1 if!l_^(10fc2 - 9fc + 2) :.'= 



fc=2 



fc(2A: + l)!!^., 



fc=3 ^ k=0 

As before, we notice that the product of the last two sums may be written as 
J^kLs o-kx'' where a^ > when k is even and a^ < when k is odd. Using this 
together with the above Taylor expansions, the fact that J^ at cos^^^{24>)d4> = 
for even k and that 

we obtain that the left-hand side of (13.2), 
32 r''^ d^ 



TlX a^((E(-l)'^'^'»'^^''^°^'(20))) cos(2</))d</> 

/•'r/2 f)2 ,,°° .. 

32/ — ^((^afc5'=cos'=+i(2(/))jjd(/.>0. 



Therefore 

3(A;' - ^;') + 2,5(3A;' + 3A;' - 2A") + 2(3A; + 3A^ - 2A') > 

and (vi) holds. 

Estimate (vii) follows now in a straightforward way: Since "iA^iS) — A{5) < 
for (5 G (0, 1/2) it is sufficient to show that 

F{S) := (1 + 2d){3Ay{5) - A{5)) - (1 - 26){3AAS) - A(6)) < 0. 

By symmetry we have F{0) = 0. Moreover, rotating n(Z, 1/2) ^ (log(2)/7r)a;z 
(Lemma 4.3) in the a;z-plane by 45°, we obtain that 3^^(1/2) —yl(l/2) = 0, proving 
(v) as well as F(l/2) = 0. By (vi) we also know that F is convex: 

F{S) < (1 - 2(5)F(0) + 2SF{l/2) = 0, 

and (vii) follows. 

Next we prove (i), that is 

d(3Ay(5) - A(5)) ^ , , , 

^ — "^^ ^-^ > for S e 1, 1/2 . 

do 
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First, 

d{3Ay{S) - A{S)) 

ds 






-'Jo OS 
As before we 


^l + 25cos{2(j)) ^(2 + (5cos(20))v/l + 2(5 
- (3 + 2,5cos(20))^/^ " (3 + 2(5cos(2<?!)))^ 

may write the integrand as 


cos(20) 


d 

85 


Vl + 2a; , ^(2 + x)Vl + 2x ^^^^ 
[ (3 + 2.)3/2+' (3 + 2x)3/2 '^"^^'^^J 


' 



d(j). 



where a; — (5cos(20). Using the Taylor series expansions calculated before it is easy 
to see that 



V2^ 

d5 

■k/2 



v/l + 2(5cos(20) ^ (2 + (5cos(20))v/l + 2,5cos(2(^) 



(3 + 2(5 cos(2<?!))) ^/^ (3 + 25 cos(20)) ^^^ 

/ IZ(-l)*''afe'^''~' cos'^'(20) + ^{-l)''+^hk6^~^ cos''+\2(p) 

•'^ '-k=0 fc=0 

where ak,bk > 0. Using J^ cos^ (2(/))d(/) = for odd j implies that 

d{3Ay{5) ^ A{5)) 
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We argue similarly to show (ii), that is 

d{3A,{5) ~ A{5)) 



> 0. 



85 



<0. 



Here 



(13.5) 



8{3A.,i5) - Ai5)) 
85 

05 

7r/2 



(1 + 4cos(20) + 2cos2(20),5)^l + 2(5cos(20) 



(3 + 2(5cos(20))3/2 



d<^ 



2cos2(2(/))(3-2(5) 



/o ^l + 2(5cos(2(^)(3 + 25 cos(2(?!)))5/2 
Substituting the Taylor expansions 

1 

Vl + 25cos(20T ^_Q 

and 



^. ^,fc(2fc-l)!! fc 



1 



(3 + 2cos(2<^))5/2 3^3^ 






,(2fc + 3)!! , 
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into (13.5) we may deduce that 
d{3A^{S) - A{5)) 
dS 
^/^ 2cos^(2(/.)(3-2(^) 



k=0 



E(-i) 



fc=0 



fe(2fc + 3)!!^fe fc,^^. 



Using that 6 < 1/2 and that all the odd terms in the product of the sums equal 
zero we see that this expression is negative, proving (ii) . 

From the fact that Ax{Q) — Ay{0) (the projection II preserves symmetry) we 
infer now that 

3A,(0) = 3/2A,(0) + 3/2AyiO) ^ A,(0) + AyiO) - / {x^ + y^)/2 

JdBin{x^+y'^>2z'^} 



< A,{0) + Ay{0) + I -z' 



z^^ -1 = A(0), 

which proves (iii). 

Combining (ii) and (iii) we obtain 3Ax{S) — A{S) < for (5 G (0, 1/2), namely 

(iv). 

Last, we verify (viii)-(ix): 



d 



85 
and 
d_ 
OS 



{3A,{d) ~ A{S)) 



(SAyiS) - A(5)) 



s=o 



.5=0 



''/^ 32 cos^ {(f>) cos(20) 16 cos(20) 
""/^ 32 sm^{(j)) cos(2(/)) 16 cos{2(j}) 



Att 



3V3 



3^3 



in 
3V3' 



D 



Proof of Lemma 4-6: 

From Theorem 4.5 we deduce that 

3A,(l/2)-A(l/2)<0, 
3Aj,(l/2)- A(l/2) = 0, and 
3A^(0) - A{0) = 3Ay{0) - A{0) < 0. 

Observe now that we have proved in (13.3) that 



3 A" + 3 A" -2 A" = -16 



It follows that 



'^/^ d^ ( ^l + 2<5cos(20) 



9^2 V 



a5n(3 + 2Jcos(2(/.))3/2 



d0>O for ^G (0,1/2). 



-{iA^{5)^A{5))-{iAy{5)-A{8))> inf [-(3A,(<5)-^(5)) -(3A,(5)-A(5))] 

(3e(0,l/2) 

=: do > for 5 £ [0, 1/2]. Now by (4.1), 

n(Zp,, 1/2) = -Ko((3A,(,5) - A{5))x^ + {iAy{5) - A{5))y^ + (3A,(5) - A{5))z-'). 

As n(Zpj , 1/2) is harmonic and thus 

{iAM - A{5)) = -((3A,(J) - A{5)) + (3A,(<5) - A{5))) > do, 
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we obtain that 

sup |Cp5 + n(Zp,, 1/2)1 > |Cp5 + n(Zp,, 1/2)1(0, 0,1) 

Bi 

= \-C-Ko{3A,{S)~A{S))\ >C + Kodo, 
and the Lemma foUows with 
(13.6) m - Kodo/2. 

Proof of Theorem J^.l: 

We begin by proving (i): For sufficiently small (5 > we have 

(13.7) 



D 



.. ^.. , , , sm'{0)smi(byded(b 

OO OO Jq Jarccot(^(l-(5)cos2(0)+5sin2(0))) 

, r^' sin^(0) sin^(0)-cos^(0) 

-'o (1 + (1- (5) cos2 (</))+ (5 sin^ (</))) ' ^'(l - (5)cos2((/)) + (5sin2(,/)) 

>_C,+4 r'^ sin^('/') sin^(0)-cos^(0) 

^^/^ (1 + (1 - (5)cos2(,/)) +,5sin2(0))^/^ ^{i _ ,5)cos2(0) +,5sin2(,/)) 

> -Ci + C2 / . ^d(j) 

I.7T/2 -, 

= -Ci+C2 



Itt/A ^COs2((/)) - (5cos(20) 

Next we notice that when (p & (37r/8, 7r/2), — cos(2(/)) > 1/V2. Consequently the 
right-hand side in estimate (13.7) is estimated from below by 

.7r/2 ^ 

Cl+C2 



S'^/s ^cos^{(t))+S/V2 



<.l-sin(37r/8) i 

> -Ci + C3 / , dt 

/-l — sin(37r/8) i i 

> -Ci + C4 / ^ ^^ - rft > -Ci + C5 log -^. 

Jo ycv CO + r Vco 

It follows that for sufficiently small 6 > 0, 



dS 
In particular. 



> -C6 + C7log(-). 



A similar calculation shows that 

ByiO)-Byi5)<(Cs~C,\og{^)y, 
so that (i) holds. 
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Next we are going to prove that 

(13.8) 
d{By{5)-B{6)) 
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■k/2 



-2cos^{(j))+6cos{2(l)) 



sin^((/)) - 008^(0) 



'0 V(l + (l-(5)cos2(0)+5sin2(0))-/ y'(l-(5)cos2(^) + <5sin2(0) 

<Cio, 

which will imply that 

(13.9) \Cy{5)~Ca{6)\<Cii 

and, when combined with (i), prove (ii) and (iv). 

In order to prove the inequality in (13.8), we make the change of variables 
cos{4>) = t, implying that 



r/2 



-2cos2(0) + <5cos(2(/)) 



sin^((/)) - 008^(0) 



V (1 + (1 - ,5) cos2(0) + Jsin2(0))^/' ^'(i _ S) cos2(0) -Ssm^{c^) 



d<i) 



-2*2 + 2,5*2 - 5 



1-2*2 



< 



(1 + t2 - 2Ji2 + J) 5/2 Vt2 - 2,5*2 + ^ ,/Y—p 



dt 



/o \/mvT~~^v^ + *2 - 25*2 

At * = 1 the singularity is of order .?"_ which is integrable. We may thus estimate 

d{By{5)-B{5)) 



dS 



< Ci, 1 



^/" J-2(l-,5)t2 
v/'5 + (1 - 2<5)t2 
1/2 -2(1 - ,5)t2 



-.dt 



-.dt 



<CiAl + V5) 

1/2 



< Ci5(l + 



) <C; 



14- 



/o V*^ + (1 - 2,5)t2 
Last, we are going to show (iii), i.e. 

B.,{5) = B.M + o{6). 

First we notice that for ,5 > 0, [Vp^l ^ on the set {ps = 0} so that we may 

differentiate Bx{S): 

(13.10) 

dBx{5) _ r'^ cos(2^)cos2(^) 



^^ "'o ^{l - 5) cos2(0) + ,5sin2(,/))(l + {l-5) cos2(0) + ,5sin2(0))5/; 

cos(20) \ / cos2(0) \ 

_(l + (l-5)cos2(^)+,5sin2(0))5/2; l,^(i_^)eos2(^)+5sin2(^)y' 



7r/2 



d0. 
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The term inside the first parenthesis is smooth for aU 6 and thus harmless. The 
term inside the second parenthesis can be estimated by 

cos^(0) ^ |cos((/))| 

^ (1-6) cosset)) + 5 sin^ict)) ~ Vl - S ' 
which is bounded for S G (0, 1/2). Using the primitive function 



C0S(2(^) COS(0) 1 ^1 - COs2((^) COS^((/)) 



C, 



(l + COs(20))5/2 ^ 2 (l+COs2(,/)))3/2 

we obtain 

B^{S)^B^{0) + SBUQ) + o{5) 

Jo (1 + cos(20))5/^ 



[lo; 

[11 

[12; 

[13 

[14 

[is; 

[16 

[17] 
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